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SUMMARY 


The equations of two-dimensional compressible flow are treated 
according to the Prandtl-Busemann small-disturbance method. In contrast 
to the usual procedure^ the independent variables are the compressible 
velocity potential and stream function and the dependent variables eire 
the rectangular Cartesian coordinates in the plane of flow. The six 
first-order differential equations corresponding to the first three iter- 
ation steps are put into complex- vector form. The particuleir integrals 
of the resulting set of three equations are then directly obtained. As 
an example, tlie general results of the analysis are applied to the case 
of subsonic compressible flow past a sinusoidal wall of small, amplitude. 


INTRODUCTION 


The problem of the integration of the equations of compressible flow 
past a prescribed solid boundary has been treated in numerous papers and 
in several diverse ways. With the exception of the hodograph- 
transformatlon method, the equations are usually solved by methods of 
successive approximations . One of these, initiated by Janzen and 
Rayleigh, starts from the incompressible-flow solution and develops the 
compressibility effects in a series of powers of the undisturbed stream 
Mach number. This method is restricted to the subsonic range, since the 
differential equations of the process are always of the elliptic type. 

In general, this method yields the best results for flows past thick 
shapes for which the critical stream Mach number is much less than unity. 
In contrast, the Prandtl-Busemann small-disturbance method of iteration 
in terms of a smn.i 1 parameter starts from the undisturbed flow, the first 
step being the well-known Prandtl linearized solution. This method is 
best suited to thin profiles for which the critical stream Mach number 
is nearly equal to \mlty and, from the beginolng, yields a good approxi- 
mation to the desired rigorous solution. Both the Janzen-Rayleigh and 
the Prandtl-Busemann methods are applicable not only to plane but also 
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to axisymmetric subsonic flovs. Tbe Prandtl-Busemann proced-ure has also 
"been utilized for the calculation of plane and axisymmetric supersonic 
flows. A good discussion of these two iterative procedures, including 
a numher of basic references, is given in reference 1. 

The author has recently considered the problem of plane inccHnpress- 
ible flow past an infinitely long, sinusoidal wall of arbitrary ampli- 
tude (ref. 2). It was found that this problem could not be treated in 
the physiceil-flow plane but had to be transferred to the plane of veloc- 
ity potential and stream fiinction. In this plane, the problem was not 
only successfully treated by the smal 1 -disturbance method but also, the 
solution was rigorously expressed in the form of a nonlinear integral 
equation. This unexpected result has prompted the present investigation 
of the form of the Prandtl-Busemann small -disturbance equations with the 
velocity potential and the stream function as independent variables. 

The equation of continuity and the condition of irrotationality form the 
basis of the procedure. Each Iteration step consists then of a pair of 
first-order differential equations of the Cauchy-Riemann type rather than 
a single second-order differential equation of the Laplace or Poisson 
type. These pairs of equations, moreover, can be put into complex- 
vector forms and their particular integrals obtained in a straightforward 
manner. Q3ie present paper contains the results of this investigation 
inclusive of the third-order step. The example of subsonic flow past a 
sinusoidal wall of small amplitude is utilized to illustrate the general 
results of the analysis. No counterpart of the integral-equation approach 
of reference 1 could be found. 


ANALYSIS 


The requli-ements of continuity and irrotationality lead to the 
relations 


where 



( 1 ) 


u,v 


velocity con^ionents in direction of x- and y-axes, 
respectively 




rectangular Cartesian coordinates in flow plane 
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p density of fluid 

density of fluid in undisturbed flow 
(f) velocity potentieil 

ijr stream fmction 


The subscripts denote differentiation with respect to the designated 
variables. The quantities x, y, (f), ijr, u, and v are nondimen- 
sional with a characteristic length associated with the solid bomdary 
as unit of length and the speed U of the undisturbed stream as unit 
of velocity. 


From Bernoulli's theorem there follows the relation 




1 



Moo^Cci^ - 1) 


7-1 


( 2 ) 


where 

M Mach number of undisturbed stream 

00 

q nondimensional speed of fluid 

7 ratio of specific heats at constant pressure and 

constant volume 


The introduction of ^ and ijr as Independent variables leads to 
the following relations in place of equations (l) : 



> ( 3 ) 



For the puirpose of deriving the iteration equations of the small- 
disturbance method, let 


X = 0 + 6Xq^ + e^X2 + e^x^ + . . . 

y = t + eyj^ + e^2 + ^^5 + . . . 
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wliere and (n = 1, 2, 3 • • • ) are functions of the Independ- 

ent variables and ijr and where e is a suitable mnn.n parameter 
which serves to regulate the course of the iteration process. 


Now, 

2 2 

I"" 2 

(x0yi[f - 


(5) 


and 


= 1 + + . . . 


= ex^^ + e^xg^ + e^XjTjf + 




y^ = 1 + eyi^ + + . . . 




( 6 ) 




Suhstituting equations (6) into eqjaatlon ( 5 ) yields 

- 1 = - 2 exi^ + €^(5x1^^ + 2 xi^yi^ + x^^^ - 2x2^) + - 

5^9j^lV^10 - + 5Xi0X20 + x^^25Z( + Wl0 + 

^1r^2t - ■ ^Alt - ^95) ■" • • • 

Also, from equation (2), 

^ = 1 - I - 1) + - 1)^ - 




48 


(8) 
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Bv means of eqmtions ( 5 ), ( 6 ), (?), ^ ( 8 )^ equations (5) ^eld the 
following sets of relations when the coefficients of the various powers 
of 6 are equated to zero: 




- yi>}r 

= -^ 152 ! 


( 9 ) 


-.2„ , _ ,r 1 M + (7 + 




= y2<ir + 1 ^ y 10 


X2\|f - “3^20 ^0^10 


( 10 ) 


>'l9P'2fl - (1 + - 

\f + (7 + 1)Mco^^19F20 + [1^^ + ■ 

|(7 + 1)M^^ + |(7 + 1) ^ j^0^ 


(11) 


= -y50 + ^ 


I yi0^ + ^10 ^10^20 






2 2 
where p = 1 - • 
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By defining 


^ + I yn = ^ 


^ + ipijr = w 


^ J yn = 


(f> - ipt = V 


and ma ki ng use of the symbolic relations 

-L = i_+ ^ 

^ ^ ^ 

^ = ipfl. . 

^ 'b^j 

equations (9), (l 0 )> and (u) can be ea^iressed in the following complex- 
vector forms: 
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Corresponding to these equations are the following second-order partial- 
differential equations: 


^iww = 0 


(15) 


_ 7 + 1 ^00 

^2ww ^ 


'2xT,,^ = - 


-■3^ 


7 + 

^ p2 


I yi/ + 


(16) 


'N 


^ + (7 + 




• 1 ^^ + ^ 10 ^ 2 ^^ 


J 


(17) 


Note that, when ^ and i}r are the independent variables, the basic 
Iteration equations (l2), (l3)^ and (l4) for x and y are of the 
first order. In contrast, when x and y are the independent variables, 
the basic iteration equations for ^ or i{f are of the second order and 
are similar to equations (l5)^ (16), and (l7)* 


CALCULATION OF PAKTICULAR INTEGRALS 


Equations (l2), (13)^ and (l^) can be Integrated without difficulty. 
The general, solutions are then obtained by the addition of arbitrary func- 
tions of w to the particular Integrals. These arbitrary functions are 
of the nature of constants of integration and are chosen in such a way 
that the bo\mdary conditions of the problem under consideration are 
satisfied. 


Both the particular and the genereil solution of equation (l2) there- 
fore is simply that zj^ is a function of w only. A particular integral 

of equation (13) can then be written as follows: 


— — 


7 + 1 ^ 
8 


|(w - 




(18) 



8 


MCA TN 5229 


By substituting this expression for Z2 into equation (l 4 ), the fol- 
lowing particular integral can he obtained; 




(w - w)z 

5 


- 2 2 

Lw ^lif + f 


2 + ^ 


+ (7 + 1) %- + 


Ij. M ^ 

J (7 + 1) %- (w - ■w)ziw^ + 3 + 3 ^ + 2(7 + 1) -2- 
0 p'^ 


i(7tl)^Vlw".2^2.p2.(7tl)^ 


i (7 + 1) ^ f ^1/^^ + 3 + + I (7 + 1) ^ 


7 + 1 ^cc 

8 r2 

I P 


f r r 

-I + Zi + (W - w)zi^J z^^ + Zi + Zi 


i (w - W)Z3_=5^ z^^l- 



NACA ™ 3229 


9 


The analysis thus fax has not heen restricted and applies equally 
well to both subsonic and supersonic flows. Actually, of course, 

and yj^ as functions of ^ and ijr are radically different for the 
two types of flow and - 1 In the case of sxipersonlc flow. In 

order to llliistrate the application of the genei^ results to a partic- 
ular flow problem, the following section contains a treatment of the 
problem of subsonic canpressible flow past a sinusoidal wall of small 
amplitude . 


SUBSONIC PTLOW PAST A SINUSOIDAL WALL 


The equation of the solid boundary Is 


y = e cos X (20) 

and the flow Is assumed to be In the direction of the positive x-axls. 
The undisturbed stream velocity U and the wave length A in I'adians 
axe utilized as units of velocity and length, respectively. The param- 
eter of the problem is e, the amplitude of the wavy wall in radians. 
The appropriate solution of equation (l2), vanishing for i}f = “, Is 
simply 




= lAe 


iw 


where A is a reail quantity to be detennined by the boun dar y condition 
at the siirface, namely, t = 0 and y = e cos x. Then 


X = ^ - eA sin ^ 

and 

cos X = pA cos 

To the zeroeth order of e, cos x » cos and, therefore, A = i, or 

X = e'^^sin 

^ p ^ 

i y = i|r + — e~P^cos ^ 

3 P P ^ 


( 21 ) 
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With Zi = ^ the appropriate general expression for Z2 can 
he obtained frcan equation (18) and is 


i 7 + ^ ^ 


4r 


8 


(Pt + l)e2i^ - 


i V 1 + 2:^ + iBe^i^ + iC 


8 p2 


4 r 2 


where B and C are real quantities to he determined hy the houndary 
condition at the surface ^ that is, ijr = 0 with y = e cos x. [Cherefore, 


■X. = (i - — sin rf - 

p y 


I 5 21^ 


sin 2J^ 


(22a) 


and 





cos X = cos <i + — < 
p 

- 

8 “ 

cos 2(f + (7 + 1) 4 


P% cos 2(f) + P^C f 


(22h) 


Frcan equation ( 22 a), to the first order of x^/p. 
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Hence, compariBon with the right-hand side of equation ( 22 b) shows that 


3% = _ 1 + 1 M ^ (v + 1 ) 


= i + y + ^ ^^£°. 

2 8 p2 


Thus, to the second order of e/p. 


I e-P'l^sin / + 1 + p2 


(p)^ ^ ^ e-2Ptsin 2/ 


(23a) 


^ y = - t + - e"P'*^cos d + l/'£'\ Ji 4. Z + 1 ?!^co „-2pilf 7 + 1 

p B 3 2lpJ 1 „2 4 :r 


, 7 + 1 i 

1 - — ^ Ptlcos 2 / 


(23b) 
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j iw / 

Again, with Z]_ = i e , equation ( 19 ) leads to the general expres- 

P 


sion for : 


z^ = 


+ i ^ ^ (1 + ^ ^ 


32 p5 




^i('w-2w) 


fl + I e-3i^ 




2 + + (y + 1 ) i (y + 1) 

p2 a p2 


M_ 


Pie 


5iw _ 


3 + P^ + 2(7 + 1) - i (7 + 1) %- 


P 


2 k 


P 


^i(2w-w) 


2 + 3^ + 


(r + 1) - i (r + 1) 


P 


2 k 


P‘ 


:.i(w-2w) _ 


1 + i P^ + 

J 


2 

9 


(7 + 1) 



-3i'vf 


+ 


4 p2 L 



where D and E are real quantities to he determined hy the boundary 
condition at the surface, that is, ijr = 0 with y = e cos x. Therefore, 



( 25 a) 
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and 


cos X = cos i — (1 

^ 2 p 


- cos 


‘ 4 ) + (1)^ - ^ M„^(l + p2)cos 39f + 


r + 1 V 

3" 


32 -2 


r- 


2 1 / 


3P^ - 3(7 + 1) ^ i (7 + 1) p 

^ pc: 


cos 0 + 


-2 + ^ 


M. 


I (7 + 1) 


cos 30 + 3^ cos 0+3^ cos 30 

J 

(25b) 


Now, from equation (25a), to the second order of e/^, 


1 6 

cos X » — — + 

2 P 


" - s(t)> ^ ^ 


cos 0 - i ^ COB 20 + 

2 p 




Comparison of this expression with the right-hand side of equation (25b) 
shows that 


b5d = - 1 - 1 p2 ^ 21+J: 

It 8 32 p2 


5 + 5P^ + 3(7 + 1) - I (7 + 1) %- 

p2 2 p2 
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P^E = 7 _ + 1 p2 _ 1 7 _ + . l ^ 1 _ 2 p2 2 ( 

24 8 24 32 p2 3 9 ^2 


Thus, to the third order of e/p. 


X = 5i( - |- e“ 3 'l^sin ^ 1 + p^ - (2p\jf + 1) e"^^'^sin 2/ + 

P 4 yp y 4 p 2 

I (7 + 1) ^ ‘ e'^^sin / - 1 p2 . i p e- 5 P^sin yj> + 

- 

2^ e“ 5 P^^ J 2 - p2 + (7 + 1) + I (y + 1) p^ sin / - 

|_ ^ P^ _ 

3 - I P^ + I (7 + 1) ^ + I (7 + 1) %- Bln 39^ - 
3 9 ^ p _ 


M ^ ^ 

4 + I (7 + 1) _S_ + I (7 + 1) ^ pt Bin 3/ 

P'^ 2 p2 


(7 + 1) ^ (pt)^8in 35^ 
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and 


7 y = i t + I e"P^cos ^ J 1 + y + ^ 

P P P 2 ^ 3 ; ] 4 p2 ® 


7 + 1 “a 


^ p" 




^ p" 


pi(f cos 20 


/ 


!■ + M 


V 


.;- i 1 r2 4 . 


7 + 1 ^ 

32 p2 




5 + 3P^ + 3(7 + 1 ) %- - i (7 + 1 ) 

2 


I- e“^'^cos 0 + 


fr + I P^']e“^^’*'cos 30 + 21±J: , 

S / 32 r2 


4 + 3P^ + 


P 


P 


3(r-n)^^-iCr-n)?^ + 5(,^i)K 

^ 2 


Pt 


p 


cos 0 - 


cos 30 + 


5 + I + I + 1) %- + I (7 + 1 ) 




Pt cos 30 + 


(7 + 1) t (p,lf)2cos 30 


( 26 t) 
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For constant values of the stream function ijr ^ 0 , equations (26) com- 
prise the parametric equations of the streamlines. The corresponding 
equations for incon^iressihle flow are obtained hy putting = 0 or 
P = 1; that is, 

X = ^ - ee“^sin 0 + i e^e'^sin 2^ + ^ e^^e~^sin 0 - e”^^sin 50)+* • • 


(27a) 

and 


y = If + ee'^^cos ^ + 


— e^(e“^cos <ji - 

8 \ 

As an example of the ccmpajrison of streamlines in incompressible 
and subsonic compressible flows, consider the case of e = 0 . 64 , 

M = 0.60 or p = 0.80, and ilf = 2 . 5 * choice of 0.64 for e would 
normally be too large if flow properties were to be calculated at the 
boundary itself with the aid of equations (26). Because of the presence 
of the exponential terms in these equations, however, a rather distant 
streamline such as i)f = 2.5 can be accxirately determined. Thus, equa- 
tions (26) become 


i “ e“^cos 20] - 

e'^’^cos 30} + . . . (27b) 


X = 0 - 0.0495 sin 0 + 0.0030 sin 20 - O.OOO8 sin 30 + • • • 


and 


y = 2.7865 + 0.0772 cos 0 - 0.0036 cos 20 + O.OOO5 cos 30 + . . . 

Equations (27), for a neighboring inccmpresslble streamline given by 
ijf = 2.6, become 


X = 0 - 0.0402 sin 0 + 0.0011 sin 20 + . . . 
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and 

y = 2.8 o48 + 0.0402 cos 0 - 0.0011 co^ 20 + . . . 

A comparison of these two sets of equations shows at a glance the well- 
loaown general result that disturbances caused by the presence of a solid 
boundary die out more slowly for compressible flow than for incompress- 
ible flow. 


CONCLUDING RBMAEKS 


Tbe use of the velocity potential 0 and the stream fimction i}f 
as Independent variables is not new in aeronautical literature, A. Olham 
and his students (partictilarly L. C. Woods) at the Oxford Engineering 
Laboratory have for years utilized Incompressible 0 and ijr as inde- 
pendent variables with the velocity vector components log i and 0 

(hodograph plane) as dependent variables. Nimierlcal relaxation methods 
were especially developed for the 0 ,ilr-plane and employed in the solu- 
tion of a mnnber of diverse problems in both Incompressible and, compres- 
sible flows, including mixed subsonic and supersonic flow. In the 
present paper, the Investigation is completely analytical and is confined 
to the physical-flow plane. In fact, the compressible velocity potential 
and stream fimction are thought of as presenting a set of orthogonal 
curvilinear coordinates associated with the shape of the solid boundary 
in the flow plane. The solution of a given flow problem is then expressed 
in the form of a pair of equations transforming the dependent rectangular 
Cartesian coordinates x and y into the Independent orthogonal curvi- 
linear coordinates 0 and ijr. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May I 7 , 195^* 
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